m ■ 
o ■ 
o 

(N ■ 
^. 
<D 

in ; 
vn ■ 


I 

o 
o 


o 
in 

cn 

O 

m 
o 


o 
o 


X 
J3 


Reply to the Comment "Intercluster correla- 
tion in Seismicity" by Helmstetter and Sornette 

The comment by Helmstetter and Sornette [1] ad- 
dresses an important point of our paper [2] that needs to 
be clarified. Any scaling measure, including also the dif- 
fusion entropy (DE) method, when applied to the earth- 
quake time series, yields anomalous scaling for two rea- 
sons, long-range time correlation (intermittency) and the 
Pareto's statistics. In fact, the DE random walker posi- 
tion is determined by the number of events (earthquakes) 
in a given time interval, which, in turns, is related to both 
the number of clusters and the cluster size. The distribu- 
tion of the time distances between clusters, ■(/'(rl™') in ref. 
[2], determines the number of clusters, while the Pareto's 
law 


(1) 


gives the probability to get a cluster of size n. The result- 
ing asymptotic scaling detected by the DE is the scaling 
of the most anomalous of the two processes, i.e. the pro- 
cess with the lowest power law exponent. 

The authors of the comment [1] discuss the Generalized 
Poisson (GP) model with an exponent 7 = 2.25 [3], which 
generates b — 0.8. In fact, according to the Continous 
Time Random Walk theory [4], the scaling is ^ = 1/(7 — 
1) if 2 < 7 < 3, (5 = 0.5 for 7 > 3 and (5 = 1 for 7 < 2. 


o GP model with Pareto exponent 7^3.5 
□ GP model with Pareto exponent 7^2.25 
A GP model with Pareto exponent y=\.lS 



FIG. 1. We plot the DE analysis for different synthetic se- 
ries generated with the GP model and three different expo- 
nents of the Pareto law, i.e. 7 = 3.5, 2.25, 1.25. The values of 
(5 obtained by fitting the asymtotic behavior are also reported 
with the corresponding error. These values increase as 7 di- 
minishes passing respectively from 0.52 ± 0.02 (for 7 = 3.5) 
to 0.78 ± 0.02 (for 7 = 2.25) and finally to 0.99 ± 0.02 (for 
7 = 1.25). 


We have repeated the calculations with 7 = 2.25 and 
the results, shown in fig.l, confirm the theoretical predic- 
tion. In the figure we plot also two other cases in order 
to illustrate how 8 changes as a function of 7. In Ref. [2] 
we explored the range 7 = 2.5 — 3. 


The experimental data yield a value 8 = 0.94 confirm- 
ing that ^/;(r['"l) has a power law tail with an exponent 
/i = 2.06. Moreover, the fact that the scaling observed 
in fig. 2 of our previous paper [2] is not influenced by 
the threshold of magnitude considered, is an additional 
evidence that what we are actually observing is the real 
scaling determined by the clusters distances. 

The LR model [2] is a minimal generalization of the 
GP model (only ■0(rt™l) changes). Therefore it is very 
interesting that a more sophisticated model, like the 
ETAS, which reproduces also other properties of seis- 
micity, yields the same 6 observed in the LR and in the 
experimental data. However we would like to point out 
that we define r^™! as the distance between two con- 
secutive clusters, while in [1] this symbol refers to dis- 
tance between two consecutive earthquakes with magni- 
tude greater than M = 5. So there is no contraddiction 
between the correlations found in [1] and our results. 

Mirko S. Mega\ Paolo AUcgrini^, Paolo Grigolini^'^''', 
Vito Latora^, Luigi Palatella^, Andrea Rapisarda^ and 
Sergio Vinciguerra^'^ 


^Dipartimento di Fisica dell'Universita di Pisa and 
INFM, via Buonarroti 2, 56127 Pisa, Italy 

^Istituto di Linguistica Computazionale del CNR, Area 
della Ricerca di Pisa, Via G. Moruzzi 1, 56124, Pisa, Italy 

^Center for Nonlinear Science, University of North 
Texas, P.O. Box 311427, Denton, Texas 76203-1427 

Istituto di Biofisica del CNR, Area della Ricerca di 
Pisa, Via Alfieri 1, San Cataldo, 56010, Ghezzano-Pisa, 
Italy 

^ Dipartimento di Fisica e Astronomia, Universita di 
Catania, and INFN sezione di Catania, Via S. Sofia 64, 
95123 Catania, Italy 

^ Osservatorio Vesuviano - INGV, Via Diocleziano 328, 
80124 Napoh, Italy 

PACS numbers: 91.30.Dk,05.45.Tp,05.40.Fb 


[1] comment by A. Helmstetter and D. Sornette, 

physics/0307134 
[2] M. S. Mega, P. AUegrini, P. Grigofini, V. Latora, L. 

Palatella, A. Rapisarda and S. Vinciguerra, Phys. Rev. 

Lett. 90, 188501 (2003). 
[3] In the literature Pareto's law is referred as the cumulative 

distribution. Our exponent 7 is related to the exponent 

reported in ref. [1] by the expression ^ — 1 + b/a. 
[4] P. Grigolini, L. Palatella, G. Raffaelli, Fractals 9, 439 

(2001). 


1 


